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I.  INTRODUCTION 


Increased  demand  for  higher  power,  light  weight,  microwave  power 
tubes  has  promoted  much  interest  in  the  design  and  development  of  more 
efficient  electron  beam  collectors.  To  this  end,  several  major  areas  of 
study  have  been  conducted  in  greater  detail.  Included  among  these  are 
conditioning  and/or  refocusing  of  the  spent  electron  beam  prior  to 
injection  into  the  collector,  multistage  depressed  collectors  (MDC),  and 
low  yield  collector  materials.  As  a  result  of  this  research,  it  is 
presently  possible  to  achieve  collector  efficiencies  upwards  of  80 
percent,  resulting  in  a  substantial  increase  in  overall  tube  efficiency. 

In  addition  to  the  areas  mentioned  above,  which  affect  changes  in 
collector  performance  via  changes  in  the  basic  collector  design,  an 
extensive  effort  has  been  made  to  develop  an  accurate  simulation  model 
of  a  circularly  symmetric  traveling-wave  tube  (TWT)  collector.  One 
model  which  has  been  used  extensively,  with  generally  favorable  results, 
was  developed  by  J.  A.  Dayton,  Jr.,  H.  G.  Kosmahl,  P.  Ramins,  and  N. 
Stanklewicz  at  the  National  Aeronautics  and  Space  Administration's 
(NASA)  Lewis  Research  Center  (LRC).  Several  papers1 ’2 ’3 ’l*  have  been 
published  verifying  the  performance  of  this  model  for  a  number  of  dif¬ 
ferent  tube  types  with  MDCs  optimized  for  a  variety  of  operating  condi¬ 
tions.  In  predicting  collector  efficiency,  the  deviation  from  actual 
results  was  less  than  10  percent  for  all  cases  reported.  In  fact,  the 
majority  of  the  computational  predictions  were  accurate  to  within  5 
percent  of  the  actual  value.  On  the  other  hand,  predictions  of 


collector  electrode  current  distributions  were  often  in  much  larger 
disagreement  with  experimentally  obtained  results. 

One  aspect  of  collector  operation  which  received  only  an  elemen¬ 
tary  treatment  in  the  LRC  model  is  that  of  secondary  electronic  emis¬ 
sion.  The  major  objective  of  this  thesis  is  to  develop  a  secondary 
emission  model  consistent  with  the  physical  emission  characteristics  of 
the  most  commonly  used  collector  materials;  an  additional  objective 
being  relatively  easy  computer  implementation. 

The  format  used  in  the  remainder  of  this  text  will  be  a  rather 
extensive  outline  and  discussion  of  the  most  pronounced  and  relevant 
emission  properties  of  metals.  This  will  be  followed  by  the  development 
of  the  emission  model.  Included  in  this  section  is  the  method  of  quan¬ 
tization  of  the  yields,  as  well  as  the  methods  used  in  the  extension  of 
previously  reported  data  to  the  entire  domain  of  primary  energy  and 
Incidence.  A  brief  discussion  of  the  method  used  in  modeling  the  TWT 
from  RF  input  to  beam  collection  will  follow.  And,  finally,  the  model's 
performance  will  be  verified  and  recommendations  for  future  areas  of 


Improvement  will  be  made 


II.  SECONDARY  EMISSION  OVERVIEW 


A.  Introduction  to  the  Overview 

Before  getting  Into  the  details  of  the  actual  emission  model.  It 
la  instructive  to  take  a  closer  look  at  some  relevant  aspects  of  second¬ 
ary  emission.  Experimental  research  into  various  modes  of  secondary 
electron  production  has  been  extensive;  however,  the  present  interest  in 
electron  collection  does  not  warrant  a  general  review  of  this  broad 
subject.  Since  many  commonly  used  collector  materials  are  polycrystal¬ 
line  in  nature,  empirical  results  for  polycrystalline  metal  surfaces, 
bombarded  by  electrons,  are  most  applicable.  Thus,  it  is  this  empirical 
data  for  which  the  following  discussion  of  secondary  emission  holds. 

The  phenomenon  of  secondary  emission  was  first  observed  by  Austin 
and  Starke  when  they  noticed  that,  given  proper  conditions,  more  elec¬ 
trons  were  emitted  from  a  sample  of  metal  than  impinged  upon  it.  This 
observation  suggested  that  some  mechanism  existed  which  promoted  elec¬ 
tron  liberation  from  the  sample.  Subsequent  experiments  have  played  a 
major  role  in  the  development  of  several  basic  concepts  presently  used 
in  the  applications  of  secondary  emission. 

B.  Secondary  Electron  Energy  Classification 

One  of  the  first  observations  made  was  that  emitted  electrons  may 
be  subgrouped  into  one  of  two  classes:  true  or  slow  secondaries,  and 
fast  or  elastically  scattered  secondaries,  with  each  group  exhibiting 
its  own  unique  emission  properties.  Figure  1  is  a  typical  curve  of 
yield  versus  secondary  electron  (SE)  energy  for  a  metal  being  bombarded 
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Fig.  1.  General  secondary  yield  for  a  metal 
as  a  function  of  secondary  energy. 

by  a  normally  incident  stream  of  primary  electrons  (PEs)  with  energy 
Ep«  Although  there  is  actually  no  clear  cut  SE  energy  which  separates 
the  elastically  from  inelastically  scattered  secondaries,  region  I  is 
generally  assumed  to  be  entirely  composed  of  slow  secondaries.  Simi¬ 
larly,  region  II  Is  assumed  to  be  entirely  composed  of  elastically 
scattered  secondaries. 

The  choice  of  the  boundary  used  to  separate  the  two  regions  is 
somewhat  arbitrary.  For  Ep  greater  than  100  eV,  some  researchers5  sug¬ 
gest  using  a  value  of  Ej,  (Fig.  1)  approximately  equal  to  1/2  Ep;  how¬ 
ever,  a  SE  energy  of  50  eV  is  more  commonly  used  as  the  transition 
point.6 ’7 *9  Therefore,  unless  otherwise  stated,  it  is  assumed  that 
any  data  discussed  in  this  thesis  were  taken  using  the  value  of  50  eV  as 
the  boundary. 
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C.  Slow  Secondary  Emission  Characteristics 


Although  treated  as  one  group,  slow  secondaries  may  be  further 
decomposed  into  two  distinct  groups  of  electrons:  true  secondaries  and 
Inelastically  reflected  primaries.  As  is  apparent  by  the  component 
names,  the  most  obvious  distinction  between  the  two  is  their  origins. 
True  secondaries  are  those  electrons  which  are  actually  liberated  by  the 
solid  being  bombarded,  whereas  inelastically  reflected  primaries  are 
primary  electrons  which  have  undergone  one  or  more  inelastic  collisions 
and  then  been  reflected  into  the  sample's  surroundings.  It  is  interest¬ 
ing  to  note  that  the  true  secondaries,  hencefore  referred  to  inter¬ 
changeably  as  slow  secondaries,  are  responsible  for  any  yields  greater 
than  unity,  since  this  is  the  only  class  of  secondaries  not  composed  of 
scattered  primaries. 

Many  experiments  have  been  performed  by  Jonker10  in  an  attempt  to 
determine  the  angular  distributions  of  three  discrete  energy  classes  of 
secondaries  emitted  by  polycrystalline  nickel.  However,  due  to  the 
range  of  primary  energies  tested,  a  fixed  threshold  of  50  eV  proved  to 
be  inadequate  in  distinguishing  the  slow  and  fast  secondaries.  Table  1 
lists  the  energy  groupings  used  by  Jonker  for  what  he  terms  the  slow, 
moderate,  and  fast  secondaries,  presented  as  a  function  of  E^.  Results 
of  these  experiments  indicate  that  the  slow  and  moderate  secondaries 
exhibit  similar  distribution  characteristics.  Thus,  it  is  assumed  that 
a  combination  of  the  two  energy  classes  produces  a  total  distribution 
which  may  be  approximated  by  the  slow  secondary  distribution. 


Table  1.  Energy  groupings  of  the  slow,  Moderate, 
and  fast  secondaries  used  by  Jonker10  in 
his  emission  distribution  measurements. 


Primary 

Electrons: 

Bombardment 

Voltage 

Slow 

Secondary 

Electrons 

Secondary 
Electrons 
with  Moderate 
Velocity 

Rapid 

Reflected 

Electrons 

25  volts 

1  to  6  volts 

10  to  15  volts 

20  to  25  volts 

100  volts 

5  to  15  volts 

45  to  55  volts 

80  to  100  volts 

450  volts 

0  to  10  volts 

* 

360  to  450  volts 

*  Current  density  was  below  observation  limit 

For  normally  incident  primaries,  with  energies  on  the  order  of  20 
eV,  Jonker  experimentally  observed  a  slow  secondary  emission  distribu¬ 
tion  very  closely  resembling  the  cosinusoidal  distribution  shown  in  Fig. 
2.  Contrarlly,  other  combinations  of  Incident  angle  and  PE  energy  were 
found  to  produce  results  which  deviated  noticeably  from  this  idealized 
distribution.  In  commenting  on  his  results,  Jonker  suggested  several 
plausible  explanations  for  this,  but  since  it  is  the  observable  changes 
in  the  distribution  and  not  the  mechanisms  giving  rise  to  those  that  are 
of  the  most  Interest  here,  these  explanations  will  not  be  given  and  the 
interested  reader  is  referred  to  reference  7.  Instead,  observations  on 
the  macroscopic  distribution  changes  with  PE  energy  and  incident  angle 


will  be  discussed 


Fig.  2.  An  Ideal  cosinusoidal  aalsslon  distribution  function  P(0). 


Fro*  Figs.  3a  and  3b,  it  is  evident  that  a  slightly  increased 
eaission,  biased  In  a  direction  approximately  antiparallel  to  that  of 
the  PE  velocity  vector,  develops  as  the  incident  angle  is  rotated 
through  45  degrees.  The  magnitude  of  emission  increase  seems  to  be  more 
pronounced  at  higher  incident  energies.  But,  until  more  detailed  dis¬ 
tribution  measurements  are  available,  no  conclusions  may  be  drawn  as  to 
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whether  or  not  additional  lncreesea  in  PC  energy  would  distort  the 
distribution  even  further.  It  way  be  inferred,  however,  that  deviations 
f row  the  cosinusoidal  distribution,  resulting  from  excitation  at  oblique 
incidence,  are  not  excessive. 

Figures  4a,  4b,  and  4c  provide  useful  insight  into  the  behavior  of 
the  slaw  secondary  emission  distribution  for  changes  in  primary  energy 
with  the  angle  of  incidence  being  held  constant.  Aesuaing  that  the  best 
approxiaation  to  a  perfect  cosinusoidal  distribution  in  each  figu.e  is 
the  25  eV  curve,  a  slight  flattening  of  the  circularly  syaaetrlcal 
distribution  occurs  with  initial  increase  in  primary  energy.  Additional 
Increases  in  energy  result  in  a  slight  elongation  of  the  initial  distri¬ 
bution  for  noraal  incidence,  indicating  a  slightly  preferred  direction 
of  ealsalon.  However,  for  the  case  of  oblique  incidence,  an  additional 
flattening  of  the  distribution  is  observed  which  indicates  e  slight 
Increase  in  the  solid  angle  of  emission  relative  to  the  surface  normal. 

Proa  Figs.  4a  and  4c,  apparently  none  of  the  energy  related  dis¬ 
tributions  deviate  froa  the  others  by  an  appreciable  amount.  Similarly, 
the  distributions  for  the  25  eV  and  the  100  eV  curves  shown  In  Fig.  4b 
also  track  each  other  reasonably  well.  Although  a  larger  deviation  from 
curves  of  slallar  angle  incidence  is  exhibited  by  the  45  eV  distribution 
of  Fig.  4b,  the  resulting  deviation,  while  significant,  is  still  not  a 
predominant  feature.  Froa  these  observations,  it  is  apparent  that 
generally  the  salssion  distributions  for  variable  energy  tend  to  track 
each  other  at  noraal  incidence  better  than  the  distributions  shown  in 
Fig.  3  for  variable  angle  of  incidence.  Thus,  it  may  be  concluded  that 
the  slow  secondary  salssion  distribution  is  slightly  more  dependent  on 


Fl«.  4.  Jonkar'a10  attaurid  angular  diatrlbut Iona  for  alow 
aacondariaa  aalttad  by  polycryatalltna  nlckal. 

(•)  t  •  0*,  dlatrlbutlona  for  aavaral  valuaa  of  E  . 

(b)  •.  •  30*,  dlatrlbutlona  for  aavaral  valuaa  of  Ep. 

(c)  Ej  ■  45*,  dlatrlbutlona  for  aavaral  valuaa  of  Ep. 
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the  eagle  of  incidence  chon  on  primary  energy  for  fine  groin  applica- 
clono. 

Overall,  Cho  moults  published  by  Jonkor  sooa  Co  lndlcoto  thot  for 
most  procdcol  applications,  Cho  slow  secondary  ealsslon  function  to  a 
first  order  approximation  Is  Independent  of  Incident  angle,  as  well  as 
prlnary  energy.  Furthemore,  this  basic  distribution  nay  be  adequately 
approximated  as  a  cosinusoidal  distribution.  Although  true  secondaries 
conprlse  s  najor  part  of  the  total  yield,  the  remaining  yield  compo¬ 
nent's  distribution  characteristics  need  also  be  examined.  In  feet,  due 
to  the  anticipated  effects  of  space-charge  suppression  on  true  secondary 
yield.  It  Is  this  class  of  secondaries  which  Is  expected  to  play  a  more 
significant  role  in  the  actual  electrode  current  distributions. 

D.  Fast  Secondary  Emission  Characteristics 

The  second  group  of  electrons  comprising  the  remaining  portion  of 
the  total  emission  Is  referred  to  collectively  as  fast,  or  elastically, 
scattered  secondaries.  This  group  of  secondaries  is  mainly  composed  of 
primary  electrons  which  have  undergone  one  or  more  elastic  collisions 
within  the  sample  end  then  been  deflected  out.  Although  only  s  smell 
fraction  of  the  total,  the  remainder  of  the  elastically  scattered  elec¬ 
trons  Is  comprised  of  lattice  electrons  which  have  undergone  an  elastic 
collision  with  a  FE  and  then  been  liberated  from  the  sample  through 
electron  substitution  or  s  similar  process. 

As  in  the  case  of  slow  secondaries,  any  differences  In  amission 
characteristics  resulting  from  the  secondaries  origin  are  Irrelevant, 


since  our  emission  model  will  treat  all  secondaries  emitted  in  the  fast 


energy  class  collectively.  Again,  the  measuresMnts  of  fast  secondary 
distributions  used  for  illustration  were  taken  by  Jonker.  In  his  report 
on  the  amission  distributions  of  polycrystalline  nickel,10  several 
•assured  fast  secondary  distributions  are  presented  for  a  variety  of  PE 
Incident  angles  and  energies*  For  convenience,  these  results  ere 
repeated  here  in  Figs.  5  and  6,  with  the  energy  range  used  for  the  fast 
secondaries,  as  before,  being  defined  in  Table  1. 

Fro*  Fig.  5,  it  is  evident  that  there  are  a  nunber  of  characteris¬ 
tics  cosmo n  to  sosm  of  these  feat  secondary  distributions  which  aerlt 
noting.  Probably  the  nost  obvious  is  the  bidirectional  nature  of  the 
distributions  shown  in  Fig.  5  for  oblique  Incidence.  For  each  of  these 
distributions,  two  distinct  lobes  are  evident,  indicating  the  existence 
of  two  local  aaission  aaxina.  Surprisingly,  emission  for  the  larger  of 
the  two  lobes  occurs  in  a  backward  direction  with  respect  to  the  compo¬ 
nent  of  the  PE  velocity  parallel  to  the  surface,  while  amission  for  Che 
smaller  lobe  occurs  in  the  forward  direction.  Thus,  simple  application 
of  Snell's  law  would  be  Inappropriate. 

Another  very  noticeable  difference  between  the  slow  and  fast 
distributions  aeasured  by  Jonker  is  the  angular  width  of  the  emission 
lobes.  As  seen  in  Figs.  S  and  6,  both  of  the  emission  lobes  for  each  of 
the  fast  secondary  distributions  are  generally  narrower  than  those  of 
the  slaw  secondary  distributions  shown  in  Figs.  3  and  4.  This  charac¬ 
teristic  indicates  that  the  fast  secondaries  till  be  emitted  over  a 
emsller  solid  angle  so  that,  statistically,  a  more  directional  emission 
will  result. 
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Fig.  5.  Jonker's10  measured  angular  distributions  for  fast 
aacondarics  omitted  by  polycryatalline  nickel. 

(a)  E„  “  100  eV,  distributions  for  several 
angles  of  incidence. 

(b)  Ep  “  450  eV,  distributions  for  several 
angles  of  Incidence. 


Fig.  6.  Jonker's1®  measured  emission  distributions  for  fast 
secondaries  emitted  by  polycrystalline  nickel. 

(a)  8^  ■  0*,  distributions  for  several  values  of  E_« 

(b)  9^  -  30*,  distributions  for  several  values  of  Ep. 

(c)  6^  ■  43*,  distributions  for  several  values  of  Ep. 
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Froa  Fig.  6,  1C  can  be  seen  that  with  Increasing  primary  energy, 
the  magnitude  of  the  secondary  lobe  decreases  and  that  of  the  primary 
lobe  increases.  Additionally,  a  decrease  in  primary  lobe  width  is  also 
observed  for  higher  energies.  The  implication  of  these  last  two  trends 
with  energy  is  that  for  a  very  high  energy  PE,  virtually  all  elastic 
emission  will  take  place  in  a  single  reverse  direction.  It  seems  fea¬ 
sible  that  this  characteristic  behavior  may  be  exploited  in  future 
collector  design,  especially  those  in  which  secondary  suppression,  not 
efficiency,  is  the  prime  concern. 

E.  Secondary  Electron  Component  Yield  Coefficients 

As  mentioned  before,  several  other  emission  characteristics  have 
been  first  observed  through  experimental  investigation.  One  of  these 
characteristics  from  which  a  complete  set  of  yield  data  may  be  obtained 
for  a  particular  material  is  the  variation  in  the  secondary  emission 
energy  distribution  with  PE  energy  and  incident  angle. 

An  indication  of  the  existence  of  the  dependency  on  primary  energy 
may  be  seen  in  Fig.  7,  which  shows  the  SE  energy  distributions  of 
molybdenum  for  a  number  of  PE  energies  impinging  at  normal  incidence. 
However,  the  utility  of  the  data  for  molybdenum  is  limited  because  the 
corresponding  data  providing  the  dependencies  upon  incident  angle  have 
not  been  reported. 

While  the  above  data  for  molybdenum  may  not  be  useful  in  itself, 
with  a  slight  explanation  it  does  imply  that  if  the  secondary  energy 
distribution  varies  with  PE  energy  and  incident  angle,  the  component 
yields  will  also.  The  integration  of  the  SE  energy  distribution  over 
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Fig. 


7.  Secondary  electron  energy  distribution  for  a  Mo  target  for 
several  primary  energies  as  taken  from  hachtnberg  and  Brauer.7 
This  figure  was  originally  presented  by  G.  A.  Harrower.11 


the  entire  SE  energy  spectrum  gives  the  total  yield,  o  ,  for  the  given 
primary  beam  energy  and  incident  angle  specific  to  that  spectrum.  Thus, 


W 

®  (E  ,  ej  -  /  i  (Ej  dE  (1) 

t  p  l  s  s  s 

o 

Similarly,  since  the  two  regions  of  Fig.  1  are  the  commonly  accepted 
domains  of  the  slow  and  fast  secondaries,  separate  integrations  per¬ 
formed  over  the  defining  energy  ranges  will  return  the  total  yield  of 
slow  and  fast  secondaries.  Thus,  letting  a be  the  total  yield  of  slow 
secondaries  and  <7^  be  the  total  yield  of  fast  secondaries,  the  relations 
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follow.  As  Is  apparent  from  defining  Eqs.  1  through  3,  the  sum  of  the 
slow  and  fast  secondary  yield  is  the  total  yield.  Thus,  the  knowledge 
of  how  any  two  of  the  above  yields  vary,  with  PE  energy  and  incidence 
for  a  given  material,  is  sufficient  in  characterizing  the  secondary 
emission  yield  properties  for  that  material. 

Using  a  spherical  retarding-f ield  analyzer,  Koshikawa  and 
Shimizu8  made  extensive  measurements  of  SE  yield  coefficients  for  smooth 
polycrystalline  copper.  In  principle,  the  analyzer  was  designed  to 
evaluate  expressions  1  and  3.  However,  due  to  the  assumed  monoenergetic 
division  of  the  slow  and  fast  secondary  domains,  slight  deviations  from 
the  physical  yield  coefficients  are  probable,  although  unimportant  for 
our  purposes. 

The  results  obtained  from  these  experiments  for  the  yield  and 
backscattering  coefficients  are  shown  here  in  Figs.  8  and  9,  respec¬ 
tively.  It  is  important  to  reiterate  that  the  data  presented  in  these 
figures  provide  complete  component  yield  information  for  polycrystalline 
copper  at  each  of  the  discrete  PE  energy-incident  angle  pairs  used  in 
the  measurements. 
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Fig.  8.  Total  secondary  electron  yield  versus  angle  of 
incidence  for  a  number  of  primary  electron 
energies  using  a  polycrystalline  copper  target.8 


Fig.  9.  Backscattering  coefficient  versus  angle  of 
incidence  for  a  number  of  primary  electron 
energies  using  a  polycrystalline  copper  target.8 


18  - 


F.  Yield  Variations  with  PE  Initial  Conditions 


There  are  several  notable  features  of  the  total  yield  curves  shown 
in  Fig.  8  which  should  be  discussed.  The  first  is  that  for  each  of  the 
primary  energy  curves  shown,  the  minimum  emission  occurs  for  normal 
incidence,  with  a  local  maximum  occurring  at  the  largest  measured 
angle.  Work  done  by  Muller12  indicates  that  the  yield  would  continue  to 
rise  for  even  larger  angles  of  incidence,  with  the  local  maximum  over 
this  range  expected  for  grazing  incidence.  From  these  results,  it 
appears  that  the  total  yield,  at  least  for  polycrystalline  copper,  is  a 
monotonically  Increasing  function  with  increasing  incident  angle. 

Another  interesting  feature  of  these  yield  curves  is  the  approxi¬ 
mately  uniform  spacing  between  curves  of  different  primary  energy.  The 
most  obvious  implication  of  this  is  that  if  the  yield  at  normal  inci¬ 
dence  is  subtracted  off,  then  the  shape  of  the  resulting  yield  curves 
with  incident  angle  will  all  be  approximately  the  same.  Thus,  it  seems 
that  this  feature  may  be  utilized  in  constructing  a  similar  family  of 
curves  for  other  materials,  once  the  variation  in  yield  with  PE  energy 
and  the  general  shape  of  the  yield  curve  versus  incident  angle  are  known 
for  these  materials. 

Additional  research*3*11*  has  also  shown  that  the  extent  of  yield 
variation  over  all  angles  of  incidence  is  very  strongly  affected  by  the 
surface  texture,  since  no  true  normal  to  the  surface  exists  for  rough 
surfaces.  Thus,  by  suitable  texturing  of  the  surface,  it  is  possible  to 
reduce  significantly  the  yield  dependence  upon  incident  angle. 
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Similar  tendencies  to  those  of  total  yields  are  exhibited  by  each 
of  the  curves  shown  in  Fig.  9  for  the  backscatterlng  coefficients.  As 
could  be  expected,  the  backscatterlng  coefficient  is  also  a  mono t on i- 
cally  increasing  function  of  the  incident  angle.  However,  a  major 
difference  between  the  two  families  of  curves  is  evident  in  the  extent 
of  functional  variation  with  Incident  angle. 

For  illustrational  purposes,  it  is  convenient  to  define  some 
figure  of  merit  for  the  total  functional  variation  with  angle.  Defining 
a  quantity  called  the  maximum  percentage  increase  (MPI)  as 

mpt  -  Functional  value  at  the  maximum  angle  .  , 
"  functional  value  at  normal  incidence^  '  ' 

it  can  be  seen  from  Table  2  that  the  MPI  of  the  total  yield  is  a  mono- 
tonically  increasing  function  of  energy,  while  that  for  the  backscatter- 
ing  coefficient  appears  to  reach  a  maximum  somewhere  between  3  and  10 
keV.  As  well  as  the  behavioral  differences  in  these  curves,  it  is 
evident  from  Table  2  that  for  all  cases  calculated,  the  MPI  of  the  total 
yield  is  greater  than  that  for  the  corresponding  MPI  of  the  backscatter- 
ing  coefficient.  This  last  observation  leads  to  several  other  interest¬ 
ing  trends  in  yield  composition. 


Table  2.  MPI  figures  for  coaponent  yields. 


0.5  keV 

1.0  keV 

2.0  keV 

3.0  keV 

6.0  keV 

10.0  keV 

Total  yield 

High 

2.16 

2.12 

1.76 

1.56 

1.25 

1.05 

Low 

1.41 

1.32 

0.991 

0.820 

0.621 

0.505 

Percent  increase 

153 

161 

178 

190 

201 

209 

Backscattered 

High 

0.577 

0.698 

0.698 

0.678 

0.633 

Low 

0.401 

0.436 

0.402 

0.380 

0.360 

Percent  increase 

144 

160 

174 

178 

176 

Percent  yield  at 
normal  inci¬ 
dence  of  slow 
secondaries 

72 

_ 

67 

51 

39 

29 

Since  It  appears  that  the  total  yield  is  greater  than  that  for 
fast  secondary  yield,  it  seems  intuitively  obvious  that  the  correspond¬ 
ing  MPI  of  the  slow  secondary  yield  should  also  exceed  chat  of  fast 
secondaries.  A  formal  proof  that  this  Is  indeed  true  is  presented  In 
Appendix  B.  By  itself,  this  observation  may  be  misleading,  since  a 
larger  MPI  does  not  necessarily  imply  a  larger  gain  in  the  number  of 
emitted  secondaries.  However,  for  the  energy  range  in  which  more  than 
50  percent  of  the  total  yield  at  normal  incidence  is  composed  of  slow 
secondaries,  changes  in  the  absolute  yield  of  slow  secondaries  are  also 
larger  than  those  of  the  fast  secondaries.  From  Figure  10,  it  can  be 
seen  that  for  polycrystalline  copper,  this  energy  range  is  approximately 
0  keV  to  3  keV. 
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A  final  observation  on  yield  versus  priaary  energy,  although  not 
very  useful  for  our  purposes,  aay  be  seen  froa  the  results  of  experl- 
aents  performed  by  Trump  and  van  da  Graaf.15  Their  results  indicate 
that  as  the  energy  of  a  normally  incident  priaary  beaa  is  increased  to 
values  on  the  order  of  200  keV,  both  the  total  yield  and  the  fast  sec¬ 
ondary  yield  asymptotically  approach  a  constant.  Since  this  trend  was 
observed  for  several  target  aetals,  it  is  reasonable  to  assume  that  this 
is  a  coaaon  characteristic  of  all  aetals. 

G.  Redesign  Universality  in  Metals 

Further  support  of  the  asymptotic  approach  of  component  yields  in 
aetals  is  the  existence  of  several  other  coaaonly  accepted  universal 
emission  properties  of  aetals.  As  lap 1 led  earlier,  included  in  this 
list  of  properties  is  the  energy  distribution  curve  shown  in  Fig.  1. 
Although  this  curve  has  several  characteristics  coaaon  to  all  aetals, 
the  aost  obvious  of  these  is  the  present  and  general  characteristics  of 
the  two  distinct  enhanced  eaission  regions. 

The  peak  located  in  the  first  of  these  regions  corresponds  to  the 
aaxlaua  yield  of  aonoenergetic  true  secondaries.  A  series  of  oeasure- 
aents  taken  over  the  energy  range  defining  region  I  of  Fig.  1  was  per¬ 
formed  by  R.  Kollath16  for  10  different  aetals.  His  results  Indicate 
that  the  difference  in  secondary  energies  at  which  this  aaxlaua  occurs 
for  aetals  is  on  the  order  of  s  2  eV,  and  approximately  centered  at  2 
eV. 

A  slightly  aore  iapresslve  result  is  the  agreesnnt  displayed  by 
aetals  over  the  entire  slow  secondary  energy  range  of  the  measured 
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energy  distributions.  Figure  11,  taken  Iron  ftachenberg  et  el.,7  repre¬ 
sents  the  spread  of  all  recorded  slow  secondary  distributions.  As  is 
evident  by  the  else  of  the  shaded  area,  a  high  degree  of  conforalty  is 
exhibited  by  all  of  the  netals  tested. 


I 


Fig.  11.  Slow  secondary  energy  spread  for  netals. 

A  second  coasson  characteristic  of  the  energy  distribution  curve  of 
Fig.  1  is  the  presence  of  fine  grain  subsidiary  naxiaa  at  slightly  lower 
energies  than  that  of  the  high  energy  naxlaua  located  in  the  second 
enhanced  emission  region.  Rudberg17  has  shown  for  several  netals  that 
the  positions  of  these  subsidiary  peaks  with  respect  to  the  peak  of  the 
high-energy  aaxlnun  are  material  specific  and  virtually  independent  of 
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primary  energy  over  the  energy  range  tested,  0  eV  <  Ep  <  400  eV. 

Che  eubeldlary  maxima  are  fine  grain  characterletlcs,  however,  their 
aentlon  here  eervee  nore  to  eapheelee  the  concept  of  metal  universality 
chan  to  provide  direct  secondary  emission  characteristics  relevant  to 
our  model. 

Perhaps  the  moat  widely-accepted  universal  emission  characteris¬ 
tics  of  metals  Is  the  energy  dependence  on  the  total  yield  for  a  nor¬ 
mally  Incident  primary  beam.  Although  the  yield  curves  for  metals 
generally  have  substantial  differences  In  absolute  terms,  when  plotted 
using  a  yield  normalised  to  the  maximum  yield  and  a  primary  energy 
normalised  to  the  primary  energy  for  which  the  stxlwm  yield  occurs, 
excellent  agreement  between  metal  results.  Figure  12,  taken  from 
Hachenberg  £t_  al. ,7  contains  measured  data  for  several  metals  after 
normalisation.  As  Is  evident  from  this  figure,  a  reasonably  tight 
apread  of  normalised  yields  is  maintained  over  the  energy  range 
plotted.  Thus,  the  normalised  yield  has  been  considered  a  universal 
property  of  metals. 

Although  extensions  of  this  property  to  primary  beaus  at  oblique 
Incidence  have  not  been  verified  experimentally  for  the  lack  of  similar 
data  corresponding  to  that  of  Fig.  12,  it  does  seem  likely  that  these 
Inferences  are  In  fact  valid.  One  line  of  reasoning  justifying  this  is 
as  follows.  As  mentioned  before,  when  subjected  to  a  primary  beam  of 
variable  Incidence,  polycrystalline  copper  experiences  an  Increase  in 
total  yield,  with  the  magnitude  of  this  Increase  dependent  on  angular 
Incidence  and  seemingly  Independent  of  primary  energy,  except  at  very 
low  energies.  Additionally,  It  should  be  noted  that  the  general  shape 
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Fif.  12.  Normalized  yield  versus  normalized  primary  energy  for  metals.  7 

of  the  yield  curve  versus  incident  angle  has  been  experimentally 
observed  to  be  similar  for  all  metals,**  and  that  the  aforementioned 
property  of  polycrystalline  copper  also  appears  to  apply  to  other  metals 
as  mail.  Since  the  major  difference  in  yield  curves  obtained  by  using  a 
normally  incident  or  an  obliquely  incident  primary  beam  is  an  approxi¬ 
mately  constant  yield  Increase  over  larger  primary  energies,  the  normal¬ 
ised  yield  curves,  although  different,  will  again  be  confined  to  toler¬ 
able  limits. 

There  are  two  other  observed  emission  characteristics  of  metals 
which  are  relevant  to  the  secondary  emission  model  being  developed.  The 
first  of  these  is  that  the  yield  of  metals  is  relatively  Independent  of 


primary  current.9  Since  thle  property  hee  been  observed  to  be  valid 
over  a  variety  of  laboratory  conditions,  no  further  discussion  of  this 
issue  is  needed. 

Unlike  the  yield  dependence  on  primary  current,  the  second  of  the 
above  referenced  characteristics  common  to  metals  has  certain  restric¬ 
tions.  Over  certain  temperature  ranges,  the  yield  of  metals  has  been 
shown  to  be  relatively  Independent  of  temperature.  Since  the  yield  is 
dependent  on  the  work  function  of  the  emitting  material,  which  in  turn 
is  related  to  surface  crystal  structure  and  impurity  concentration,  this 
property  does  not  hold  over  temperature  changes  for  which  either  of  the 
two  preceding  parameters  has  been  significantly  altered.  These  may  not 
be  the  only  restrictlona  in  order  for  this  property  to  be  valid,  but  at 
present  the  author  is  unaware  of  any  others. 

Although  this  discussion  has  by  no  means  covered  all  aspects  of 
secondary  emission,  it  has  covered  most  concepts  needed  to  develop  a 
reasonable  model.  With  this  in  mind,  attention  will  be  focused  on  the 
task  of  developing  a  physically  accurate  secondary  emission  model. 
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III.  SECONDARY  EMISSION  MODEL 


A.  General  Considerations  of  the  Emission  Model 

Fro*  Che  preceding  discussion  of  the  secondary  emission  character¬ 
istics  of  metals  due  to  electron  bombardment,  there  are  certain  obvious 
emission  properties  which  must  be  considered  in  the  development  of  any 
emission  model.  The  major  of  these  being: 

1.  Yield  and  component  yield  variations  with  primary  energy  and 
incident  angle. 

2.  Emission  distribution  variations  with  primary  energy  and 
Incident  angle. 

3.  Variation  of  the  modeling  process  (due  to  both  controllable 
and  uncontrollable  operating  conditions). 

Although  all  of  the  preceding  properties  of  secondary  emission  are 
important  considerations  in  the  model  being  developed,  probably  the  most 
Important  consideration  is  the  range  of  PE  conditions  for  which  the 
model  is  expected  to  perform.  This  is  true,  since  this  range  not  only 
plays  a  major  role  in  the  model's  overall  range  of  validity,  but  also 
the  choice  of  modeling  process  used  for  the  emission  properties  previ¬ 
ously  mentioned. 

Since  the  emission  model  will  be  used  primarily  in  TWT  simula¬ 
tions,  it  seems  reasonable  to  require  that  the  model  be  able  to  accu¬ 
rately  accommodate  incident  angles  in  Che  range  0°  <  0^  <  90* ,  as  well 
as  primary  energies  in  the  range  of  0  keV  <  Ep  <  10  keV.  The  upper 
limit  of  10  keV  was  chosen,  since  few  TWTs  operate  with  a  cathode-to- 
helix  voltage  exceeding  10  keV,  and  those  that  do  generally  operate 
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their  collectors  at  depressions  such  that  the  beam  is  depressed  below 
the  10  keV  threshold  when  collected.  In  addition,  the  basic  secondary 
emission  data  are  not  available  above  10  keV.  Although  the  model  needs 
to  be  the  most  accurate  in  the  domain  defined  above,  provisions  should 
also  be  made  for  PEs  having  initial  conditions  outside  of  this  domain. 

Acknowledging  the  trade  off  in  model  generality,  much  of  the 
secondary  emission  model  will  be  based  on  empirical  data  in  order  to 
accommodate  the  above  range  of  PE  Initial  conditions.  Since  these  data 
are  material  dependent,  a  choice  of  target  material  is  necessary. 
Because  copper  is  one  of  the  most  commonly  used  TWT  collector  materials 
in  industry  today,  empirical  data  for  polycrystalline  copper  will  be 
used  when  available.  Instances  in  which  appropriate  data  are  not  avail¬ 
able  for  a  copper  target,  similar  data  for  a  substitute  material  will  be 
used.  Fortunately,  this  is  not  expected  to  present  a  problem,  since 
many  referenced  sources  do  use  copper  as  the  emitting  material. 

B.  Component  Yield  Model  Development 

Two  secondary  emission  parameters,  which  have  been  observed  to 
change  drastically  over  the  range  of  PE  Incident  angles,  are  the  total 
and  backscattered  yields.  As  mentioned  before,  both  are  monotonically 
increasing  functions  with  total  yield  increases  on  the  order  of  90 
percent  being  typical  for  a  60°  change  in  incident  angle,  and  similar 
results  applying  to  backscattered  yields.  Since  the  use  of  textured 
surfaces  in  order  to  reduce  the  component  yield  dependencies  on  incident 
angle  is  uncommon  in  today's  TWT  industry,  the  variation  of  yields  with 


-  29  - 


f 


incident  angle,  as  mentioned  before,  should  be  included  in  the  emission 
model. 


Strong  dependencies  have  also  been  observed  for  the  yield  compo¬ 
nents  versus  primary  energy.  Published  results8  show  that  for  poly¬ 
crystalline  copper,  the  total  yield  may  almost  triple  over  the  energy 
range  0.5  keV  <  Ep  <  10  keV,  with  the  backscattered  coefficient  decreas¬ 
ing  nearly  50  percent  over  that  same  range.  As  in  the  case  of  yield 
cofficients  with  primary  incidence,  this  relation  should  also  be 
Included  in  the  model  being  developed. 

The  curves  of  Figs.  8  and  9,  taken  from  Koshikawa  et  al.,8  serve 
as  the  primary  source  for  all  data  necessary  to  develop  the  required 
functional  relationships  existing  between  the  component  yields  and  PE 
initial  conditions  for  polycrystalline  copper.  Rather  than  developing  a 
two-variable  expression  relating  these  quantities,  a  number  of  single 
variable  expressions  were  developed  for  fixed  values  of  the  second 
variable.  Since  the  developed  expressions  are  valid  over  a  continuous 
domain  of  the  first  variable  only,  a  further  interpolation  over  the 
second  variable  is  generally  necessary  for  any  given  set  of  PE  initial 
conditions. 

Since  the  yield  curves  appear  to  be  better  behaved  with  angle  than 
with  energy,  the  yield  coefficients  were  represented  by  using  the 
primary  energy  as  the  continuous  variable,  and  the  incident  angle  as  the 
interpolating  variable.  Additionally,  since  PE  initial  conditions  are 
only  approximate,  as  well  as  the  fact  that  the  yield  curves  are  gener¬ 
ally  smooth  over  incident  angle,  any  higher  order  interpolation  over 
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incident  angle  would  provide  no  better  approximation  than  a  linear 
interpolation. 

A  possible  further  benefit  of  using  the  linear  interpolation, 
although  not  originally  intended  as  such,  deals  with  the  sign  of  the 
predicted  error.  Since  all  yield  curves  shown  exhibit  the  same  nearly 
parabolic  form,  calculated  values  using  linear  interpolation  will  be 
slightly  higher  than  values  taken  directly  from  the  curves.  Combining 
this  observation  with  the  fact  that  Kowhikawa  and  Shimizu  mentioned  that 
their  measured  backscattered  coefficients  for  large  angles  of  incidence 
may  have  been  lower  than  actual  yields  due  to  the  possible  presence  of 
stray  fields  in  the  measurement  apparatus,  it  can  be  seen  that  the 
linear  interpolation  may  result  in  a  slightly  more  accurate  backscat¬ 
tered  coefficient. 

In  the  case  of  the  total  yield,  however,  the  increase  in  predicted 
value  related  to  the  Interpolation  Is  not  expected  to  compensate  for 
experimental  errors,  since  these  measurements  are  considered  by  the 
experimenters  to  be  accurate.  Thus,  this  process  would  return  a 
slightly  more  erroneous  slow  secondary  coefficient.  However,  due  to  the 
relatively  low  energy  levels  of  the  true  secondaries  and  the  relative 
strength  of  the  retarding  fields  in  most  collection  regions,  the  error 
Introduced  in  the  slow  secondary  yield  is  expected  to  be  less  detri¬ 
mental  to  model  operation  in  contrast  to  beneficial  compensation  noted 
for  the  fast  secondary  yield. 

The  angular  step  size  chosen  between  successive  component  yield/ 
primary  energy  expressions  was  5  degrees.  The  primary  motivation  behind 
this  choice  of  step  size  was  that  it  seemed  to  offer  a  reasonable 


trade-off  between  compactness  and  interpolation  accuracy,  assuming  that 
the  above  assumption  does  not  apply  and  that  the  data  presented  in  Figs. 
8  and  9  are  valid. 

As  is  evident  from  the  preceding  discussion,  the  structure  of  the 
above  component  yield  model  should  allow  for  reasonable  results,  since 
the  error  introduced  by  this  method  of  data  storage  and  retrieval  should 
be  relatively  minor  compared  to  the  error  introduced  by  other  aspects  of 
collector  simulation  modeling  to  be  discussed  shortly.  Although  the 
structure  of  the  model  is  sound,  several  conditions  must  be  met  in  order 
to  assure  reasonable  performance  when  this  structure  is  implemented. 
The  first  of  these  is  that  the  data  published  by  Koshikawa  and  Shimizu 
are  accurate,  and  the  second  is  that  a  reasonable  fit  to  the  data  may  be 
made  over  both  PE  energy  and  incident  angle. 

Although  the  accuracy  of  Figs.  8  and  9  are  slightly  beyond  our 
control,  a  rough  idea  as  to  the  reliability  of  the  model  due  to  the  data 
used  may  be  obtained  by  comparison  with  other  experimental  results. 
Using  the  results  of  several  experiments,  McKay9  indicated  that  most 
likely,  maximum  yield  of  copper  for  a  normally  incident  beam  is  approxi¬ 
mately  1.3  occurring  for  a  PE  energy  of  about  0.6  keV.  The  0.5  keV 
yield  curve  presented  by  Koshikawa  et  al.  suggests  that  the  total  yield 
at  normal  incidence  is  approximately  1.4,  with  the  corresponding  maximum 
yield  possibly  higher  still.  Since  there  are  a  number  of  subtle  mate¬ 
rial  variables  which  may  directly  affect  the  predicted  yields,  while  not 
making  the  measurements  valid,  this  10  percent  discrepancy  in  the  maxi¬ 
mum  yields  is  not  unrealistic  or  unacceptable  in  terras  of  the  overall 
collector  model. 
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Although  this  characteristic  is  not  intended  as  conclusive  proof 
that  the  data  shown  in  Figs.  8  and  9  are  accurate,  this  along  with 
several  other  predicted  trends,  with  both  incident  angle  and  energy, 
does  indicate  that  these  data  are  at  least  realistic.  These  trends  will 
become  more  obvious  once  the  component  yield  functional  relationships 
are  developed.  However,  in  lieu  of  further  evidence  either  proving  or 
disproving  the  validity  of  Figs.  8  and  9,  the  data  provided  by  Koshikawa 
et  al.  will  be  assumed  accurate  in  the  remainder  of  the  thesis. 

A  second  condition,  which  must  be  met  in  order  for  the  component 
yield  model  to  operate  reasonably,  is  an  accurate  extrapolation  of  the 
reported  data  to  the  domains  outside  of  those  measured.  Since  measured 
data  were  only  reported  for  incident  angles  in  the  range  0°  <  0^  <  65°, 
with  the  larger  angles  of  incidence  changing  most  rapidly,  the  region 
65°  <  0^  <  90°  may  contain  significant  extrapolation  error.  However, 
since  all  yield  curves  presented  seem  to  be  smooth  and  follow  an  almost 
quadratic  relation  over  the  incident  angle  range  measured,  a  relation 
which  from  other  researchers  is  expected  to  hold  over  larger  angles  of 
incidence  as  well,  appreciable  error  is  not  expected  to  be  much  of  a 
problem  in  the  extrapolated  region. 

For  each  yield  curve  shown  in  Figs.  8  and  9,  yield-angle  data 
pairs  were  supplied  to  a  least  squares  program.  The  order  of  the  fitted 
equations  was  clamped  to  be  at  least  one  lower  than  the  number  of  data 
pairs  being  fitted;  and  two  types  of  equations,  a  polynominal  in  X  and 
an  exponential  in  X,  were  used  as  the  form  of  the  approximating  equa¬ 
tion.  Invariably,  as  expected,  all  predicted  best-fit  equations  were 
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polynomials  in  X,  with  generally  only  minor  contributions  made  by  third 
and  higher  order  terms. 

The  resulting  yield  relations  plotted  over  all  incidence  are  shown 
in  Figs.  13  and  14.  The  dashed  curve  in  Fig.  13  was  taken  from  refer¬ 
ence  12  and  represents  a  total  yield  curve  for  incident  angles  approxi¬ 
mately  reaching  78  degrees.  Although  the  PE  energy  used  in  generating 
the  curve  is  unknown,  the  similarity  in  shape  to  the  adjacent  curves 
serves  to  Indicate  that  the  method  used  for  extrapolating  the  yield  at 
larger  angles  of  incidence  gives  reasonable  results. 

To  complete  the  component  yield  model,  equations  relating  the 
component  yields  to  the  primary  energy  were  calculated  in  much  the  same 
way  as  those  used  for  extrapolation  to  large  incident  angles.  Two 
additional  features  were  employed  in  the  least  squares  fit  routine  in 
order  to  yield  physical  curve  shapes.  The  first  of  these  was  to  fit  the 
data  to  a  polynomial  in  1/X,  a  relation  which  seemed  applicable  in  view 
of  the  previously  mentioned  asymptotic  character  that  physical  yield 
curves  exhibit  at  larger  primary  energies.  The  second  of  these  was  the 
method  used  in  choosing  the  best  fit  equations. 

A  smoothing  routine  based  on  the  null  hypothesis18  was  also  added 
to  the  least  squares  fit  routine  and  a  new  set  of  equations  was  calcu¬ 
lated.  The  resulting  set  of  best  fit  equations  was  all  polynomials  in 
1/X,  found  to  be  superior  to  polynomials  in  X.  Several  of  the  curves 
generated  from  these  quations  are  shown  in  Figs.  15  and  16. 


Fig.  13.  Total  yield  curves  extrapolated  to  grazing  incidence  using 
the  results  of  Koshikava  and  Shimizu8  and  compared  to  data 


found  in  reference  12  and  originally  reported  by  Muller 


12 
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Fig.  14,  Backacattered  yield  curvet  extrapolated  to  grazing 

incidence  ualng  the  reaulta  of  Koahlkava  and  Shimizu.0 
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Fig.  16.  Backscattered  yield  curves  generated  using  the  backscattered  yield  expressions  developed  for 
the  emission  model.  These  expressions  were  generated  using  the  data  presented  In  Fig.  14. 


Certain  features  of  these  approximating  curves  appear  to  make  this 
family  of  component  yield  curves  a  better  choice  for  the  component  yield 
model.  The  first  of  these  is  the  general  shape  of  the  total  yield 
curves  shorn.  As  mentioned  in  the  previous  discussion,  not  only  is  this 
a  universal  characteristic  of  metals,  but  the  relative  shape  is  virtu¬ 
ally  independent  of  incident  angle.  A  second  characteristic  of  the 
approximating  family  of  curves  which  merits  noting  is  the  smoothness 
over  the  energy  range  plotted. 

C.  Component  Yield  Model 

Although  the  preceding  characteristics  would  seemingly  indicate 
that  representing  the  component  yields  by  polynomials  in  1/X  provides  a 
good  representation  of  physical  values,  characteristics  of  these  curves 
also  exist  which  limit  the  range  of  their  validity,  of  which  the  most 
obvious  occurs  for  small  values  of  primary  energy. 

Because  the  approximating  functions  are  polynomials  in  1/X,  com¬ 
ponent  yields  approach  infinite  values  as  the  primary  energy  approaches 
zero.  Since  the  yield  must  be  finite,  this  representation  is  obviously 
invalid  for  very  low  primary  energies,  and  thus  a  special  adaptation  is 
used  in  this  range. 

A  second  range  of  primary  energies  for  which  the  derived  relations 
are  invalid  is  for  high  primary  energies.  This  is  due  to  slight  dis¬ 
crepancies  in  asymptotic  values  of  total  yield  and  backscattered  yield 
for  the  incidence  of  5 degrees  which,  if  left  as  is,  would  predict  a 
small  negative,  slow  secondary  yield.  Thus,  this  requires  that 
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Modification  also  be  made  In  the  coaponent  yield  aodel  for  high  primary 
energies  as  well. 

Deterioration  in  the  curve  fits  to  the  data  was  observed  below 
0.75  keV.  To  be  conservative,  the  derived  equations  were  assumed  to  be 
valid  over  the  range  0.75  keV  <  Ep  <  10  keV.  Note  that  the  upper  bound¬ 
ary  of  10  keV  was  not  established  to  be  the  upper  Halt  for  a  consistent 
set  of  yield  equations,  but  rather  assumed,  since  empirical  data  were 
available  only  to  that  point.  The  resulting  energy  ranges  with  the 
component  yield  functions  used  in  each  region  are  shown  in  Fig.  17. 

Several  points  can  be  aade  about  the  component  yield  model. 
First,  the  difference  between  the  upper  and  lower  boundaries  of  the 
cross-hatched  region  at  any  given  primary  energy  is  the  slow  secondary 
yield  for  a  given  primary  energy.  Expressed  in  a  different  way, 

o  *  o  +o  (5) 
t  8  t 

This  point  should  be  obvious  from  Eqs.  1-3. 

A  second  point  of  Interest  is  the  alternate  method  used  to  deter¬ 
mine  the  yield  coefficients,  and  thus  the  component  yields  for  primaries 
with  energies  less  than  0.75  keV.  As  is  obvious  from  the  figure,  the 
fast  secondary  yield  is  assumed  to  be  constant  for  this  region,  assuming 
the  same  value  at  0.75  keV.  In  this  same  region,  the  total  yield  tapers 
linearly  from  the  backscattered  value  at  a  primary  energy  of  0  keV  to 
the  predicted  value  at  0.75  keV. 


|1.0  2.0  3.0  0.0  5.0  6.0  7.0  8.0  9.0  1040  11.0 

0.75  Ep(keV) 


Fig.  17.  Tbs  three  distinct  regions  of  yield  modeling  and 

their  respective  values  for  total  and  backscattered 
yields  as  used  In  the  secondary  emission  model. 
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Several  things  activated  this  choice  of  modeling  for  this 
region.  First,  while  a  finite  yield  is  known  to  exist  in  the  limit  as 
the  primary  energy  approaches  0  keV,  the  limiting  value  is  unknown. 
Additionally,  it  is  also  known  that  the  yield  for  very  low  energies  is 
composed  of  a  high  percentage  of  fast  secondaries,  with  this  percent¬ 
age  increasing  even  further  for  further  decreases  in  primary  energy. 
This  last  trend  continues  until  a  primary  energy  of  approximately  10  eV 
is  reached,  for  which  case  the  total  yield  and  the  fast  secondary  yield 
are  approximately  the  same.6  As  can  be  seen  from  this,  the  model  and 
the  actual  phenomenon  are  qualitatively,  if  not  quantitatively,  in 
agreement  for  this  region. 

A  final  point  about  the  component  yield  model  is  the  yield  treat¬ 
ment  for  primary  energies  above  10  keV.  For  primary  energies  in  this 
range,  the  component  yields  are  assumed  to  have  reached  their  asymptotic 
value,  and  are  thus  treated  as  constants.  Even  though  the  physical 
asymptotes  are  not  actually  approached  until  much  larger  primary  ener¬ 
gies,  the  approximation,  as  is  the  case  of  the  low  primary  energy  model¬ 
ing,  is  not  expected  to  increase  the  error  substantially  in  the  overall 
collector  model,  due  to  the  relatively  small  number  of  electrons  being 
collected  with  energies  in  this  range  for  most  TWT  collector  designs. 
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D.  Component  Yield  Quantization 


Along  with  the  component  yields,  the  emission  model  should  include 
an  emission  angle  dependence.  This  should  be  obvious  in  light  of  the 
previous  discussion  on  emission  distributions.  As  before,  much  of  this 
modeling  will  be  based  on  empirical  data  available  for  polycrystalline 
nickel. 

Since  the  model  is  to  be  included  in  a  computer  simulation,  some 
discrete  representation  of  the  observed  continuous  distribution  is 
necessary.  As  mentioned  earlier,  the  observed  distributions  of  fast  and 
slow  secondary  yield  vary  significantly,  especially  with  solid  angle  of 
emission.  It  vas  thus  decided  that  the  broader  of  the  two  distribu¬ 
tions,  that  of  the  slow  secondaries,  was  to  be  represented  by  three 
rays,  and  that  the  relatively  narrow  fast  secondary  distribution  would 
be  represented  by  a  single  ray.  The  limit  of  four  secondaries  per 
primary  was  a  prudent  one  in  view  of  the  enormous  memory  requirements 
and  increased  execution  time  for  each  additional  ray. 

As  mentioned  before,  results  published  by  Jonker  show  that  it  is 
reasonable  to  assume  a  cosinusoidal  distribution  for  true  secondaries, 
as  shown  in  Fig.  2.  Since  three  rays  are  to  represent  the  true  second¬ 
ary  emission,  it  Is  workable  to  ascribe  to  each  ray  one-third  of  the 
true  secondary  emission  current.  However,  since  the  program  used  in  the 
actual  collector  simulation  is  a  two-dimensional  code,  the  division  of 
the  total  current  into  thirds  will  be  restricted  to  the  projected  planar 
current  and  not  the  actual  three-dimensional  current.  Since  the 
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collector  electrodes  are  assumed  to  be  axially  symmetric,  this  approach 
should  provide  a  good  first  approximation  to  the  actual  distribution. 

Appendix  A  contains  a  detailed  derivation  of  the  expression  for 
the  total  planar  current,  with  emission  angles  between  0^  and  6^*  Since 
a  closed-form  solution  could  not  be  found,  all  integrations  needed  to  be 
performed  numerically.  The  numerical  quadrature  method  chosen  for  all 
integrations  was  the  trapezoidal  rule  with  Richardson  extrapolation. 

The  first  integration  performed  was  done  over  all  emission  angles 
so  that  the  total  planar  current  yield  was  known.  This  value  was  then 
used  as  a  normalization  factor  for  subsequent  integration. 

Since  the  division  of  current  elements  into  thirds  provides  a 
method  for  distributing  the  secondary  yield,  but  not  the  emission  angle 
of  each  component,  it  was  decided  that  each  third  would  be  emitted  at 
the  median  angle  of  its  respective  wedge  section.  By  definition,  one  of 
these  would  be  the  direction  normal  to  the  surface,  however,  the  other 
two  are  not  so  obvious.  Thus,  an  iteration  employing  bisection  of  the 
Integration  range  was  used  in  determining  these  remaining  launch 
angles.  The  resultant  discrete  distribution  of  slow  secondaries  used  in 
the  emission  models  is  shown  in  Fig.  18. 
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Fig.  18.  Launch  angles  for  the  quantized  yields  used  In  the  secondary  emission  model. 

Launch  angles  1,  2,  and  3  represent  the  three  emission  angles  for  slow  second¬ 
aries.  Launch  angle  4  represents  the  emission  angle  for  the  backscattered  yield 


Determining  the  emission  angle  for  the  fast  secondary  yield  is 
much  more  straightforward.  Table  3  contains  a  tabulated  version  of 
Jonker’s  measurements  for  fast  secondary  distributions.  In  many  of  the 
distributions  presented  by  Jonker,  the  distribution  maximum  serves  as  an 
approximate  axis  of  symmetry  as  well;  thus,  the  angle  of  maximum  emis¬ 
sion  was  chosen  as  the  single  direction  of  fast  secondary  emission.  For 
most  of  the  values  presented  in  the  table,  the  emission  angle  is  approx¬ 
imately  2/3  of  the  incident  angle,  regardless  of  primary  energy.  Thus, 
a  phenomenological  adaptation  is  made  for  setting  the  fast  secondary 
launch  angle,  despite  the  lack  of  adequate  physical  explanation  for  its 
occurrence. 


Table  3.  Tabulated  fast  secondary  emission  angle  as  taken  from  Jonker. 
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E.  Emission  Model  Stability  with  Operating  Conditions 

Since  TWTs,  as  well  as  a  number  of  other  electron  devices  which 
require  electron  collection,  are  subject  to  a  variety  of  operating 
conditions,  it  is  important  to  examine  the  effects  of  the  emission 
model's  flexibility  and  reliability  with  changes  in  operating  condi¬ 
tions.  For  several  reasons,  primary  current  and  temperature  may  vary 
considerably  in  both  time  and  space  for  any  given  device  or  collector. 
Since  these  parameters  are  not  easily  controllable  by  appropriate  col¬ 
lector  design,  the  effects  of  these  on  the  model  are  relevant,  and  thus 
their  impact  as  related  to  TWT  collectors  will  be  discussed. 

By  their  very  nature,  TWTs  modulate  the  current  and  velocity  of  an 
electron  beam  in  extracting  energy  from  it.  Noting  that  the  beam  cur¬ 
rent  varies  for  different  applications  as  well,  a  further  scaling  of 
primary  current  entering  the  collector  may  result.  Since  the  current 
entering  the  collector  is  time  dependent,  it  is  intuitively  obvious  that 
collected  primary  currents  would  also  differ  significantly  with  both 
time  and  space.  As  discussed  before,  however,  metals  exhibit  a  univer¬ 
sal  tendency  to  maintain  relatively  constant  yield  characteristics  with 
PE  current,  thus  variaitons  in  PE  current  are  not  believed  to  be 
directly  responsible  for  errors  in  the  emission  model. 

What  is  true,  however,  is  that  temporal  and  spatial  variations  in 
primary  current,  as  well  as  the  ambient  temperature  of  the  device,  set 
up  temperature  gradients  over  the  collection  surface.  Pulsed  devices 
working  on  the  principle  of  electron  beam  switching  further  complicate 
matters  in  that  they  may  make  these  gradients  a  stronger  function  with 
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time.  As  before,  however,  it  has  been  shown  that  yield  characteristics 
of  metals  are  virtually  independent  of  temperatures.  All  of  this  leads 
to  the  conclusion  that  temperature  and  current  densities,  which  are  not 
easily  controllable  in  practice,  fortunately  are  not  expected  to 
adversely  affect  the  model's  performance  in  simulating  physical  collec¬ 
tors. 

Although  other  operating  conditions  may  exist,  which  do  in  fact 
affect  the  performance  of  the  emission  model,  these  are  not  apparent  to 
the  author.  Thus,  it  is  concluded  that  no  discrepancy  between  actual 
and  modeled  yield  data  will  be  solely  the  result  of  operating  condi¬ 
tions.  The  final  consideration  of  the  secondary  emission  model  to  be 
discussed  is  its  range  of  validity. 

F .  Range  of  Validity 

Since  much  of  the  emission  model  previously  developed  has  its 
basis  in  empirical  data,  certain  conditions  must  be  satisfied  in  order 
for  the  emission  model  to  be  valid.  The  first  of  these  is  that  whatever 
the  target  material  being  used,  it  must  exhibit  similar  yield  tendencies 
versus  incidence  and  primary  energy.  The  second  is  that  the  surface  be 
of  comparable  smoothness  as  the  sample  used  in  obtaining  the  data. 
Finally,  the  target  material  should  be  a  metal  (or  else  operating  condi¬ 
tions  may  change  the  performance  of  the  model). 

If  the  above  conditions  are  met,  a  further  embellishment,  a  yield 
scaling  factor,  has  been  added  to  the  emission  model.  For  similar 
metals  (i.e.,  copper  and  nickel),  this  feature  appears  to  be  useful; 
however,  extension  to  other  suitable  collector  electrode  materials  may 
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be  Inappropriate  more  often  than  not.  For  these  metals,  a  more  elabo¬ 
rate  normalization  of  the  curves  would  prove  worthwhile.  Possible  areas 
of  further  adaptation  of  this  model  may  involve  an  inclusion  of  normal¬ 
ized  yield,  as  previously  discussed.  Alternately,  a  normalization  of 
yield  data  with  angle  of  incidence  may  be  explored,  possibly  by  using 
the  MPI  figure  previously  defined,  or  some  other  method.  These,  how¬ 
ever,  will  not  be  Included  in  our  emission  model;  and  thus,  the  validity 
of  our  model  is  limited  to  the  restrictions  previously  mentioned. 


IV.  COLLECTOR  ANALYSIS  MODEL 


A.  Introduction  to  the  Analysis  Model 

As  nay  be  seen  from  the  previous  discussions,  an  improved  second¬ 
ary  emission  model  by  itself  is  not  enough  to  ensure  good  MDC  simula¬ 
tion.  Not  only  is  this  true,  but  the  addition  of  a  more  accurate  emis¬ 
sion  model  may  actually  make  analytical  analysis  of  collectors  less 
accurate,  resulting  from  an  increased  significance  of  PE  incidence  and 
energy.  This  point  demonstrates  the  need  for  certain  inclusions  in  the 
total  collector  analysis  model  to  ensure  at  least  a  reasonable  represen¬ 
tation  of  PE  trajectories  upon  collection. 

Due  to  the  generally  favorable  success  of  the  LRC  model,  the 
overall  approach  of  collector  simulation  used  in  that  model  will  also  be 
used  in  this  one.  As  a  consequence  of  this,  the  simulation  will  be 
approached  as  a  two-stage  process;  the  first  of  these  being  the  modeling 
of  the  beam  from  the  RF  input  of  the  TWT  to  the  entrance  of  the  collec¬ 
tor,  or  the  Interaction  region,  and  the  second  being  the  simulation  of 
the  beam  in  the  MDC  region. 

For  the  most  part,  the  modeling  of  the  individual  regions  will 
also  be  paralleled  with  that  of  the  LRC  model.  Thus,  with  only  few 
exceptions,  all  features  of  the  LRC  models  are  also  Included  in  the 
models  used  in  this  simulation.  Additionally,  several  other  features, 
not  specifically  acknowledged  as  being  included  in  the  LRC  model,  have 
also  been  incorporated  into  the  RF  Interaction  simulation. 
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B.  RF  Interaction  Region  Modeling 


For  Che  purposes  of  determining  a  PE  initial  condition  profile 
upon  entrance  into  the  MDC,  a  large  signal  program  which  computes  both 
RF  power  and  beam  data  was  employed .  A  version  of  the  large  signal 
program  employed  by  Dayton  jst_  al.  was  used  for  this  RF  region  model¬ 
ing.  In  way  of  a  brief  description,  this  simulated  beam-RF  interaction 
was  based  on  J.  E.  Rowe's*9*2®  nonlinear  interaction  theory  developed 
for  a  traveling -wave  amplifier.  Other  considerations  provided  for  the 
LRC  model  are: 

1.  The  introduction  of  RF  circuit  losses. 

2.  Provisions  for  attenuators  and  severs. 

3.  User  definable  phase  velocity  and  circuit  impedance. 

4.  User  definable  variable  circuit  radius. 

5.  Simple  and  complex  PPM  focusing. 

6.  Relativistic  correction  to  the  axial  equation  of  motion. 

7.  The  use  of  the  complete  second  order  circuit  equations  for 
the  RF  wave. 

8.  The  retention  of  all  terms  involving  C  (the  square  of 
Pierce's  gain  parameter). 

9.  Optional  number  of  deformable  disks  or  rings  used  in  modeling 
the  electron  beam. 

10.  The  use  of  the  hot  phase  velocity  for  the  beam  loaded  RF 
wave. 

11.  The  optional  inclusion  of  a  backward  wave. 
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In  addition  to  the  above,  the  large  signal  program  used  in  this 
modeling  also  provided  for: 

12 .  The  optional  inclusion  of  multiple  forward  waves. 

13.  The  relativistic  correction  to  all  equations  of  motion. 

14.  An  energy  conservation  algorithm  to  facilitate  an  energy 
conserving  electron  beam. 

One  of  the  major  influences  on  the  trajectories  of  the  individual 
electrons  comprising  the  electron  beam  is  the  total  power  extracted  from 
it.  Thus,  the  difference  between  measured  extracted  power  and  that 
predicted  by  the  large  signal  gain  program  was  used  as  a  measure  of  how 
well  the  predicted  spent  beam  approximated  the  physical  one.  An  inter¬ 
action  of  inputs  to  the  large  signal  gain  program  was  performed  until  an 
accurate  extracted  power  was  achieved.  Throughout  the  course  of  this 
interaction,  only  parameters  assumed  not  to  be  known  accurately  were 
varied.  Included  in  these  are  the  phase  velocity,  attenuation,  and 
interaction  impedance. 


Collector  Region  Modeling 


Simulation  of  the  collector  region  was  accomplished  by  using  a 

slightly  modified  version  of  an  electron  gun  program  originally  devel- 

21 

oped  by  N.  J.  Dionne  and  J.  J.  Krahn  for  the  Naval  Research  Labora¬ 
tory.  In  turn,  this  program  had  its  origins  in  the  Kirstein-Hornsby 
electron  gun  program. 

The  approach  used  to  generate  electron  trajectories  is  to  approxi¬ 
mate  the  path  traveled  by  the  charge  groups  as  a  series  of  straight  line 
segments  in  two  dimensions.  This  ballistic  approach,  which  involves  the 
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determination  of  the  kinematic  equations,  is  then  used  in  determining 
the  charge  group's  new  position  assuming  a  finite  travel  time.  The 
process  is  then  repeated  until  the  charge  group  is  collected. 

The  major  features  of  this  collector  model  include: 

1.  The  use  of  maganetostatic,  electrostatic,  and  space-charge 
force  terms  (assumed  axially  symmetric)  in  the  calculation  of 
all  kinematic  equations  of  motion. 

2.  A  finite  travel  time  step,  small  in  relation  to  the  RF  period. 

3.  The  use  of  the  secondary  emission  model  developed  earlier. 

4.  An  energy  conservative  treatment  of  the  collection  of  pri¬ 
maries. 

5.  An  option  to  run  fast  secondaries  only,  slow  secondaries  only, 
both  species  of  secondaries,  or  no  secondaries. 

6.  The  optional  inclusion  of  multiple  bounces  of  secondaries  when 
any  single  species  option  of  secondaries  is  employed  (slow 
secondaries  are  then  represented  as  a  single  current  ray 
Instead  of  the  three  current  rays  used  when  both  species  of 
secondaries  are  employed). 

Incorporation  of  the  space-charge  forces  into  the  collector  model 
raises  several  questions  as  to  the  validity  of  an  underlying  assumption 
implicit  in  the  collector  code.  Since  the  beam  in  the  gun  region  is  not 
bunched  axially,  it  is  reasonable  to  assume  that  a  dc  current  exists  for 
gun  analysis.  Thus,  it  was  assumed  that  the  discrete  charge  groups  may 
be  replaced  by  continuous  current  rays.  Although  reasonable  for  gun 
analysis,  this  approximation  does  not  strictly  hold  for  collector  simu¬ 
lation. 
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Since  Che  space-charge  forces  may  contain  some  RF  component  due  to 
the  velocity  modulated  beam,  the  use  of  a  continuous  current  ray  will  at 
best  approximate  the  time  averaged  RF  space-charge  effects  on  predicted 
trajectories.  However,  it  is  felt  that  the  use  of  this  resultant  quasi¬ 
steady-state  solution  will  provide  a  sufficiently  accurate  simulation. 

One  final  point  of  interest  is  the  conversion  of  the  time  depen¬ 
dent  output  of  the  large  signal  gain  program  to  the  steady-state  condi¬ 
tions  required  for  the  collector  simulation.  All  time  dependence  was 
eliminated  by  an  axial  shift  of  all  charge  rings  to  a  common  plane. 
During  this  shift,  however,  the  magnetic  flux  referenced  back  to  the 
cathode  is  not  necessarily  conserved.  Thus,  appropriate  measures  were 
also  taken  in  calculating  the  radial  coordinates  to  ensure  that  the 
resulting  data  still  conserve  magnetic  flux. 
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V.  COLLECTOR  MODEL  PERFORMANCE 


A.  Single  Bounce  Secondary  Model  Performance  with  Varying  Spike 

The  accuracy  of  the  collector  simulation  model  was  examined  by 
comparison  of  predicted  and  measured  collected  current  distributions  for 
an  existing  continuous  wave  TWT.  The  TWT  used  for  verification  incor¬ 
porated  a  two-stage  MDC,  with  all  collecting  surfaces  being  either 
nickel  or  nickel  plated .  One  concern  in  using  this  tube  was  that  its 
geometry  ensured  that  a  slightly  nonsyrametric  magnetic  field  existed  in 
the  collector  region. 

Other  basic  operating  characteristics  of  the  TWT  are  as  follows: 


Operating  bandwidth 
Beam  current 
Beam  voltage 
Maximum  output  power 
Beam  focusing 
First  stage  depression 
Second  stage  depression 
Spike  depression 


7.6  GHz 
0.240  A 
10.15  kV 
275  W 
PPM 

49  percent 
72  percent 
100  percent 


All  verification  runs  were  made  at  the  single  frequency  of  10.4 
CHz.  Additionally,  a  saturated  mode  of  tube  operation  was  employed  in 
all  of  the  simulations. 

The  angular  and  kinetic  energy  distributions  of  the  modeled  satu¬ 
rated  beasi  used  are  shown  in  Fig.  19.  Evident  from  the  figure  is  that  a 
reasonably  large  range  of  particle  trajectories  exists  at  this  entrance 
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KINETIC  ENERGY  (VOLTS) 


(a)  Energy  distributions  given  as  a  percentage 
of  the  original  beam  kinetic  energy. 

SPENT  BEAM  ANGULAR  CLASS  DISTRIBUTION 
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(b)  Angular  distributions  given  as  the  angle  between 
the  electron  trajectory  and  the  axis  of  the  TWT. 

Fig.  19.  Spent  beam  distributions  upon  entrance  into  the  collector 
simulation  program  used  for  model  verification. 
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plane.  Consequently,  an  even  wider  range  of  trajectories  is  expected  to 
exist  at  the  points  of  collection.  By  utilizing  this  enlarged  domain  of 
the  secondary  emission  model,  a  reduction  in  the  cumulative  error  of 
collected  current  distributions  introduced  by  the  emission  model  for  any 
individual  ray  should  be  realized. 

Although  in  the  preceding  context,  statistical  averaging  offers 
more  insight  into  overall  model  performance,  limitations  to  this 
approach  also  exist.  The  most  notable  of  these  is  the  strong  RF  compo¬ 
nent  of  beam  current  also  associated  with  a  saturated  mode  of  operation 
and  neglected  in  our  modeling  process. 

For  a  10. A  GHz  signal,  the  periodicity  of  the  beam  current  is 
approximately 


Tfe  *  1/10.4  *  109  *  9.615  x  10  11  s 

Also  noting  that  the  axial  dimension,  da>  Qf  the  collector  being  modeled 
is  on  the  order  of  0.06  meters,  the  transit  time  for  a  typical  electron 
in  the  collector  may  be  roughly  approximated  as 

T  -  *  d  /v  *  d  //TnV  3  1  x  10-9  s 
trans  a  avg  a 

By  virtue  of  the  relative  magnitudes  of  these  times,  a  significant 
number  of  RF  beam  bunches  may  be  present  in  the  collector  at  any 
instant.  Thus,  predicted  current  distributions  may  be  considerably 
different  from  experimentally  observed  distributions. 
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To  mitigate  the  consequences  of  this,  verification  was  performed 
using  curves  rather  than  a  point  by  point  comparison.  The  initial 
verification  was  generated  by  varying  the  spike  potential.  Since 
changes  in  the  spike  potential  do  not  greatly  alter  the  fields  where  the 
bunching  effect  is  most  pronounced,  the  entrance  to  the  collector,  the 
error  in  computed  trajectories  due  to  the  omission  of  RF  beam  bunching 
in  our  modeling,  should  remain  somewhat  constant  over  the  range  of  spike 
potentials  plotted.  Correlations  between  measured  and  predicted  curve 
shapes  should  then  provide  a  measure  of  how  well  the  emission  model 
performs. 

Figure  20  shows  the  current  distributions  obtained  by  three 
methods;  one  method  being  measured  data  and  the  other  two  constructed 
from  simulations. 

In  order  to  make  the  measurements  on  an  unpackaged  tube,  the  tube 
was  operated  at  a  greatly  reduced  duty  cycle  of  0.1  percent.  A  10  us 
pulse  width  was  also  used  to  ensure  that  any  transient  beam  effects  are 
negligible  in  the  measurements. 

The  initial  runs  incorporated  a  yield  scaling  factor  of  1,  since 
nickel  and  copper  have  similar  maximum  yields  at  normal  incidence.9  In 
each  of  the  curves  generated  for  the  individual  stages,  some  qualitative 
similarity  to  measured  results  is  evident  above  spike  voltages  of  2500 
volts • 

Although,  in  this  range,  the  magnitudes  of  this  set  of  generated 
curves  differ  dramatically  from  the  empirical  curves,  both  sets  of 
curves  exhibit  maxima  and  minima  at  approximately  the  same  locations. 
The  most  interesting  of  these  being  the  minimum,  which  occurs  in  the 
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(b)  First  stake  current  versus  spike  depression 


Pig.  20.  Continued 


SPIKE  current  vs.  spike  depression 

cb-ieise  v  Ebi-siaa  v  Eba-ae-ia  v 


(d)  Spike  current  versus  spike  depression 


spike  current  lg.  For  this  region,  both  experimental  and  predicted 
results  Indicate  that  the  spike  acts  as  an  emitter. 

For  the  range  of  spike  voltages  below  2500  Vt  however,  predicted 
results  do  not  appear  to  agree  quite  as  well  with  measured  distribu¬ 
tions,  either  quantitatively  or  qualitatively.  This  is  most  evident  in 
the  body  and  first  stage  currents,  IW8  and  Ibl»  respectively,  where 
current  variations  in  excess  of  45  percent  of  an  expected  constant  value 
are  predicted. 

Consideration  of  the  magnitude  of  the  spike  current,  Ig,  in  the 
region  where  the  spike  acts  as  an  emitter  suggests  that  the  simulations 
may  be  overoptlmistic  in  their  predicted  yields.  A  second  set  of  veri¬ 
fication  curves  was  generated  using  a  scaling  coefficient  of  0.7  in 
order  to  examine  this  possibility.  The  results  of  these  simulations  are 
represented  by  the  remaining  set  of  curves  shown  in  Fig.  20. 

As  in  the  first  set  of  distributions,  there  are  obvious  similari¬ 
ties  between  the  shapes  of  the  predicted  and  measured  curves,  especially 
for  lower  values  of  spike  depression.  Additionally,  although  not 
readily  observable  from  the  curves  presented  in  this  form,  similarity  to 
measured  distributions  hold  better  at  higher  spike  depressions  as  well. 

Summing  the  body  and  first  stage  currents  yields  the  curves  shown 
in  Fig.  21.  As  in  che  case  of  the  individual  currents,  measurements 
indicate  a  reasonably  constant  current  for  spike  depressions  greater 
than  80  percent.  On  the  other  hand,  predictions  made  using  a  scaling 
factor  of  1  indicate  that  a  significant  variation  in  the  combined  cur¬ 
rents  over  this  range  of  spike  depression  will  result,  this  variation 
being  on  the  order  of  20  percent  of  the  measured  value.  However,  for 
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Fig.  21.  Combi 


similar  curve*  generated  using  a  scaling  factor  of  0.7 ,  this  vsrlstlon 
is  reduced  to  roughly  10  percent.  Furthermore,  the  absolute  accuracy  of 
the  prediction*  in  this  range  was  also  significantly  improved. 

There  are  several  possibilities  why  a  reduced  yield  proved  to 
yield  slightly  better  results  than  the  yield  which  would  seemingly  be 
accurate  for  the  nickel  collection  surfaces  of  the  collector  used,  the 
most  probable  being  the  macroscopic  texture  of  the  collecting  surface. 

Current  and  Fox13  have  shown  that  through  suitable  surface  textur¬ 
ing,  the  yield  of  the  emitting  material  may  be  reduced.  The  target 
surface  used  in  obtaining  the  copper  yield  data  was  first  polished 
mechanically,  electrolytical ly ,  and  finally  cleaned  by  argon-ion  bom¬ 
bardment,  resulting  in  a  surface  texture  similar  to  that  shown  in  Fig. 
22a.  On  the  other  hand,  the  machined  parts  making  up  the  collecting 
surfaces  would  probably  have  a  surface  texture  similar  to  that  shown  in 
Fig.  22b. 

A  second  possibility  is  that  the  secondary  lobe  previously  dis¬ 
cussed  as  being  observed,  but  ignored  in  the  backscattered  distribution, 
does  in  fact  govern  an  appreciable  amount  of  the  backscattered  cur¬ 
rent.  Since  the  net  effect  of  including  this  forward-directed  current 
would  almost  certainly  be  a  net  decrease  in.  the  current  deposited  at 
lesser  depressed  stages,  this  is  also  a  very  plausible  explanation. 
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(a)  Sputtered  copper  surface  (b)  As-machined  copper  surface.  X1000. 
having  cones.  X1000. 


(c)  Sputtered  copper  surface  having  needle-like  protuberances. 

Fig*  22.  Scanning  electron  micrographs  of  copper  surfaces  showing  texture. 


And  finally,  a  third  source  of  error  aay  be  Introduced  by  the 
emisalon  of  higher  order  secondary  yield.  Figure  23  shows  the  predicted 
trajectories  for  each  pair  of  the  points  comprising  the  0,7  yield  scal¬ 
ing  distributions.  It  should  be  evident  from  these  figures  that  at 
least  some  portion  of  the  tertiary  and  higher  order  secondary  yield  will 
propagate  toward  higher  depressed  stages,  effectively  reducing  the  yield 
if  only  secondaries  are  considered  in  the  analysis,  at  least  for  larger 
values  of  spike  depression. 

Note  that  the  oalsslon  of  the  higher  order  secondaries  may  also  be 
the  primary  cause  of  the  body  and  first  stage  currents  for  depressed 
spike  potentials,  since  for  this  case  there  is  less  of  a  retarding  field 
in  the  vicinity  of  secondary  collection,  thus  allowing  backscattered  as 
well  as  true  tertiary  and  higher  order  secondaries  to  be  emitted.  Since 
the  alow  secondaries  are  low  energy  electrons,  it  is  almost  assured  that 
they  will  either  be  recollected  somewhere  on  the  emitting  electrode  or 
at  even  less  depressed  stages.  This,  coupled  with  the  relative  size  of 
the  slow  secondary  yield,  could  lead  to  a  significant  secondary  migra¬ 
tion  to  less  depressed  stages. 
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(a) 


Eb  -  10,150  V,  Em  -  5180  V,  Eb2  -  2840  V, 
and  Es  “  0  V,  all  with  respect  to  cathode* 
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(b)  Eb  -  10,150  V,  EbJ,  -  5180  V,  Efc2  -  2840  V, 
and  Eg  “  1030  V,  all  with  respect  to  cathode* 

23.  Simulated  paths  of  primaries  and  secondaries  for  a  yield 
scaling  coefficient  of  0.7  and  all  electrode  potential 
combinations  used  in  generating  Figs.  20  and  21. 
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(c)  Eb  -  10,150  V,  Ebl  -  5180  V,  Efe2  -  2840  V, 
and  Eg  ■  2000  V,  all  with  respect  to  cathode* 


(d)  Eb  -  10,  150  V,  Eb,  -  5180  V,  Eb,  -  2840  V, 
and  Eg  “  2600  V,  all  with  respect  to  cathode* 


Fi«*  23*  Continued 
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(f)  Eb  .  10,150  V,  Ebi  -  5180  V,  Eb2  -  2840  V, 

•nd  Eg  "  2790  V,  all  with  respect  to  cathode. 

Fig.  23.  Continued 
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(g)  Eb  -  10,150  V,  Ebl  -  5180  V,  Eb2  -  2840  V, 

and  E6  ■  2840  V,  all  with  respect  to  cathode. 

Fig.  23.  Continued 
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B.  Single  Bounce  Secondary  Model  Performance  with  Varying  Second  Stage 

A  second  set  of  curves  was  generated  by  holding  the  spike  at  90 
percent  depression  and  varying  the  second  stage  depression  from  50  to  70 
percent.  For  continuity,  both  yield  scaling  factors  of  0.7  and  1  were 
used.  The  resultant  current  distributions  are  shown  in  Fig.  24. 

While  both  sets  of  curves  exhibit  a  noticeable  improvement  in 
qualitative  agreement  to  measured  results  over  the  entire  energy  range 
plotted,  the  most  notable  occurs  for  a  yield  scaling  coefficient  of  0.7 
and  when  the  body  and  first  stage  currents  are  treated  collectively. 
For  this  case,  the  qualitative  agreement  between  measured  and  predicted 
curves  is  very  good  for  all  cases  except  the  spike. 

While  the  error  in  the  spike,  body,  and  first  stage  currents  is 
probably  due  to  the  omission  of  higher  order  secondary  yields  in  our 
modeling,  the  higher  overall  current  in  the  lesser  depressed  stages  is 
most  likely  due  to  too  large  a  yield  scaling  factor.  This  may  be  due  to 
either  of  the  two  reasons  previously  mentioned. 
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AXES  ORIGIN 


Fig.  24.  Men 


current  versus  second  stage  depression 


(c)  First  stage  current  versus  second  stage  depression. 


(d)  Second  stage  current  versus  second  stage  depression. 


(e)  Spike  current  versus  second  stage  depression. 


VI.  CONCLUSIONS 


Overall,  Che  collector  simulation  suffered  from  the  same  major 
inaccuracy  as  the  LRC  simulations,  namely,  incorrect  predictions  of 
lesser  depressed  stages.  The  fact  that  the  predictions  exhibited  some 
of  the  characteristic  shapes  of  the  measured  curves  would  seemingly 
indicate  that  the  model  followed  changes  in  primary  initial  conditions 
reasonably  well. 

The  sources  of  error  discussed  earlier  in  the  text  are  probably 
the  areas  in  most  need  of  improvement.  At  the  time  of  this  writing,  the 
Raytheon  Company  is  presently  working  on  a  model  which  employs  a  3-ray 
representation  for  both  the  true  as  well  as  backscattered  yield.  Pre¬ 
liminary  results  are  in  better  agreement  to  actual  results  than  the  1- 
ray  backscattered  model  employed  in  this  thesis. 

The  problem  of  incorporating  all  higher  order  secondaries  is 
workable,  however,  at  great  expense  with  today's  computer  technology. 
For  situations  in  which  only  1  species  of  secondaries  is  emitted,  the 
model  developed  in  this  thesis  is  adequate.  Out  of  all  of  the  curves 
generated  in  this  thesis,  this  condition  would  only  apply  to  the  first 
point  of  each  of  the  curves  in  Fig.  20.  Since  these  points  happen  to  be 
in  excellent  agreement  to  the  measured  results  anyway,  there  appears  to 
be  very  little  benefit  in  employing  this  option.  The  curves  for  which 
both  fast  and  slow  secondaries  were  emitted  yielded  much  more  appreci¬ 
able  error  than  secondaries  only,  and  thus  were  omitted  from  the  text  of 
this  thesis. 


» 
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The  final  problem  of  surface  texturing  may  be  addressed  by  suit¬ 
able  application  of  the  scaling  factor.  This,  however,  does  not  account 
for  any  preferential  emission  directions  which  may  also  be  a  consequence 
of  the  surface  texture.  Again,  while  this  is  being  written,  the 
Raytheon  Company  is  presently  undertaking  the  inclusion  of  an  angular 
offset  which  may  prove  useful  in  optimizing  model  performance  given  a 
fixed  material  and  surface  texture.  Until  more  work  is  reported  in  this 
area,  however,  these  two  features  prove  little  more  than  an  additional 
degree  of  freedom  until  experience  dictates  otherwise. 

Finally,  other  areas  of  research  which  may  be  considered  as  well 
are  the  extension  of  these  yield  curves  and  distributions  to  any 
metal.  Several  possible  methods  of  doing  this  have  already  been  dis¬ 
cussed. 

Although  much  work  still  needs  to  be  done  on  the  model  to  achieve 
consistent  quantitative  agreement  for  all  cases,  the  expense  is  "but  a 
drop  in  the  proverbial  collector  bucket"  compared  to  savings  in  future 
MDC  design. 
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APPENDIX  A 


DERIVATION  OF  PLANAR  CURRENT  EXPRESSION 

Since  Che  purpose  of  this  derivation  is  to  eventually  have  an 

expression  for  the  planar  current  in  the  y-z  plane,  it  is  necessary  to 

♦ 

find  expressions  relating  the  initial  current  weighting  vector,  Vj ,  to 
the  y  and  z  component  vectors,  V2  and  V^,  respectively.  These  vectors 
are  shown  in  Fig.  A.l  for  an  arbitrary  direction  defined  by  the  angles  6 
and  ♦ . 

♦ 

The  vector  Vj  may  be  broken  up  into  components  as 

*1  “  *3  +  ^4  "  l*3l  *  +  $4  (A*1) 


Vector  $ 

4 


may  also  be  written  as  the  sum  of  its  component  vectors. 


♦ 


^2 + 


(A. 2) 


Also,  V^|  and  |  may  be  related  to  | V ^ |  by  the  relations 


\\\  -  |vj  sin(0) 


(A. 3) 


and 


V3|  -  jvj  cos(0) 


(A. 4) 


» 
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Similarly,  che  magnitudes  of  ^  and  ^  may  be  related  to  |  as 


l^2l  “  1^4 1  8in^) 


(A. 5) 


lv5l  *  |vj  C08(*) 


(A. 6) 


Substituting  Eq.  A. 4  into  Eq.  A.l  gives 


V  -  |V  |  cos(0)  z  + 


(A. 7) 


Now,  using  Eqs.  A. 2,  A. 3,  A. 5,  and  A. 6,  Eq.  A. 7  may  be  written  as 


^  *  |$j|  cos(0)  z  +  |^j|  sin( 0 )  sin($)  y  +  |$j|  sin(0)  cos($)  x  (A. 8) 


Breaking  $  up  into  components,  with  one  component  directed  in 

A 

the  x  direction  and  one  component  normal  to  that  direction  (y-z  plane), 
we  may  write 


V1  "  Vlx  x  +  Vli  1 


(A. 9) 


Equation  A. 18  may  now  be  put  into  the  form  of  Eq.  A. 9  by  defining 


vu  5  |VJ  (cos2(0)  +  sin2(0)  sin2(*)^ 


(A. 10) 


and 


1  =  jvj  (sin(0)  sin($)  y  +  cos(0)  z)/Vu  (A. 11) 

By  definition,  |vj  |  is  just  the  distribution  function,  which  describes 
the  number  of  electrons  emitted  per  unit  time  per  solid  angle. 

As  mentioned  previously,  this  distribution  function  takes  the  form 
of 

jVj|  *  Jo  cos(0)  (A. 12) 

From  Eqs.  A. 10  and  A. 12,  we  see  that  the  magnitude  of  the  total  number 
of  emitted  electrons  per  unit  area  per  unit  time  projected  on  the  y-z 
plane  is  simply 

1  2  2  2  \1/2 
-  Jo  cos(0)  fcosv0)  +  sin  (0)  sin  (#)J  (A. 13) 

Thus,  the  total  number  of  electrons  emitted  per  unit  time  in  the  y-z 

plane  may  be  determined  by  multiplying  V  by  the  solid  angle  and  inte¬ 
grating  over  a  specified  region  defined  by  0  and  The  unit  solid 

angle  may  now  be  written  as 

d?  -  sin(0)  d@  d*  (A. 14) 

An  expression  for  the  differential  planar  current  emitted  for  a  given 
direction  may  now  be  written  as 
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1/2 


dJ  -  Jo  cos(0)  (cos^(0)  +  sin^(0)  sin^($)^  sln(0)  d0  d4>  (A. 15) 

y-z  \ 


Due  to  symmetry,  integration  need  only  be  performed  over  the  hemisphere 
defined  by  x-z  plane.  An  expression  for  the  total  planar  current  con¬ 
tained  between  the  two  angles,  0^  and  0^,  may  now  be  derived  from  Eq. 
A. 15  as 

9  2  *  /  ,  0  0  X1/2 

/  dJ  -Jo  /  cos(0)  sin(0)  /  (cosz(0)  +  sin^C©)  sinz($)J  d<>  d9 

y~Z  9 i  0  K 

(A. 16) 


since 


2  2 
sin  (x)  -  1  -  cos  (x) 


(A. 17) 


Eq.  A. 16  may  be  written  as 


\l/2 


/  dJy_2  “Jo  /  cos(0)  sin(0)  /  ^cos^(0)  +  sin^(0)(l  -  cos^(4*)]j  d<p  d0 

(A. 18) 


which  may  also  be  written  as 


z  "  /  o  o  \l/2 

/  dJ  “Jo  /  cos(0)  sin(0)  /  fl  -  sin^(0)  cosz($)J  d$  d0  (A. 19) 
y~z  a  n  '  ' 


Defining 
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4  =  »/2  -  v 


(A. 20) 


we  see  thac 

d$  •  -dv  (A. 21) 

and  noting  that  the  limits  also  change,  Eq.  A. 19  may  now  be  written  as 

9  2  -ir/2  \l/2 

/  dJ  -  -Jo  /  cos(0)  sin(0 )  /  (l  "  sin2(0)  cos2(if/2  -  v)J  dv  d0 
yz  e  ir/2  V 

(A. 22) 

Since 

cos(ir/2  -  x)  -  sin(x)  (A. 23) 


Eq.  A. 22  becomes 


0 


fd J  -  J1  /  cos(0) 

y-2  0L 


x/2  f  i  \  ^ 2 

sin(0)  /  (1  -  sin  (0)  sin  ( y ) 1  dv  d0 

-tt/2V 


The  integral  on  v  may  be  evaluated  by  noting  that 


n 


(» 


dv 


E( v ,  a) 


Equation  A. 24  may  be  written  as 


-  87  - 


* 


MICROCOPY  RESOLuflON  ^EST  ^CHART  . 
- BUREAU  Of  STWjOAROS^l^-A 


'2  |V  -  */2 

/  dJ  -Jo  /  cos(0)  sin(8)  E(v,  sln(0))  (A. 26) 

1  y**  a  v  '  lv  -  -w/2 


v  -  -t/2 


Evaluating  the  term  inside  the  integral,  Eq.  A. 26  becomes 


/  dJ  -  Jo  /  cos(6)  sin(0)  {E(v/2,  sin(0))  -  E(-ir/2,  sin(0))}  d0 
y'2  \ 


(A. 27) 


It  can  also  be  shown  that  E(-v,  k)  -  -E(v,  k),  so  Eq.  A. 27  may  be 
written  as 


/  dJ  -  Jo  /  cos(0)  sin(8)  {e(*/2 ,  sin(0))  +  E(*/2,  sin(6))} 
^  8i 


2Jo  /  cos(0)  sin(0)  E(it/2,  sin(0))  d0 

9i 


(A. 28) 


Noting  that  E(ir/2,  k)  is  also  a  special  case  of  E(v ,  k) ,  we  will  define 
the  special  case  as  E(k)  and  also  note  the  relation  that 


/o  !  „2,  2  l2  *  3  .4  (2n  -  1)M  k2n  _ 

E(k)  —  v/2  1  •  2  k  —  .  s  k  ~  ...  —  ■  (on  —  ••• 

22  •  4 2  2n  n!  (2n  1} 

(A. 29) 


Using  this  fact,  Eq.  A. 28  may  be  rewritten  as 


/  dJ  -  vJo  /  cos(0)  sin(0) 


,  _  y  f(2n  -  mij2 

,n  .  (2n 

n-1  2  nl 

„  0 


(2n  -  1) 


(A. 30) 


Since  Che  integral  of  Che  sum  is  the  sum  of  the  integrals,  Eq.  A. 30  may 
be  rewritten  as 


/  dJ  ■  xJo  </  cos(d)  sin(O)  d8 


V  11 2  ~ 

" .  n 

n-1  2  n ! 


<2n  -  1) 


/  cos(0)  sin2n+1(8)ld8  (A.31) 


In  general. 


/  cos(0)  sin(0)  d0  ■  — sinB+i(0) 

IS  T  1 


(A.32) 


So,  substituting  Eq.  A.32  into  Eq.  A.31,  we  get  an  expression  for  the 
projected  planar  current  emitted  from  the  region  defined  by  0^  and  e2  as 


J  -  x Jo  J 1/2  sin2(e2)  -  sin2^) 


<2^1  (sin2(n+1>(02)  -sin2(n+1)(01))[(A.33) 


n*»l|  2  nl 


Equation  A. 33  may  now  be  solved  using  numerical  procedures. 
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APPENDIX  B 


PROOF  OF  RATES  OF  INCREASE  WITH  INCIDENT  ANGLE 


Definitions 


i(y  »i) 

*(V  ei] 

6i} 

stEp,  et) 


is  the  current  of  the  primary  beam  with  energy  Ep  and 
making  an  angle  of  0^  with  respect  to  the  surface 
normal 

is  the  total  secondary  yield  for  the  primary  current 

l(Ep,  e1) 

is  the  backscattered  yield  for  the  primary  current 

i(Ep,  »*) 

is  the  true  secondary  yield  for  the  primary  current 

Ky  y 


y(Ep,  0^)  is  the  total  yield  coefficient  for  the  primary  current 

>(y  y 

b(Ep,  Oj  is  the  backscattering  coefficient  for  the  primary 
current  l(Ep,  Bj) 

with  9 ^  shown  in  Fig.  B.l  as  being  the  angle  of  incidence  with  respect 
to  the  surface  normal. 

For  brevity,  define  a  general  function  of  the  form  f(E  ,  0  )  to  be 

b  i 

f_  0..  Now  for  normal  Incidence  by  definition, 


®E  0  "  XE  0  bE  0 
P  P  P 


(B.l) 


also 
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Fig.  B.l.  Definition  of  incident  angle  as  used  in  this  appendix 
and  everywhere  in  the  rest  of  this  report. 


XE  0  yE  0 
P  P 


(B.2) 


and 


*E  0  fyE  0  "  bE  0^ 
P  V  P  P  / 


(B.3) 


for  a  beam  with  incidence  of  67  degrees,  we  have 


®E  67  *  1E  67  bE  67 
P  P  P 


also, 


(B.4) 


YE  67  "  XE  67  bE  67  * 
P  P  P 


(B.5) 


and 
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SE  67 (yE  67  "  bE  67) 
P  \  P  P  7 


(B.6) 


Defining  Che  percenc  Increase  of  a  functional  value  with  angle  as  the 
ratio  of  the  functional  value  at  the  largest  angle  to  the  functional 
value  at  normal  incidence,  we  may  write  the  percent  Increase  in  true 
secondary  yield,  fS,  as 


SE  67  *E 


F  2 


67  (yEp67  ~  bEp67 

v  cv  ■  v) 


) 


(B.7) 


P  ~P~ 

and  the  percent  lncreaae  in  backscattered  yield,  Fb,  as 


BE„67  IE  67  bE  67 

JE _ E _ 


B 


Eo° 


I„  .  b  . 
E  0  eO 
P  P 


(B.8) 


Now,  in  order  to  ensure  that  the  only  change  in  yields  is  that  due  to 
the  angular  dependence,  let 


IE„67  ~  h  0 
P  P 


(B.9) 


With  the  above  relation,  Eqs«  b«7  and  B>8  become,  respectively. 


„s  (V7'V7) 

(EP°  EP°j 


(B.10) 


and 


P 


B 
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Defining 


rE_67 


rE_0 


we  a«y  rewrite  Eq.  B.10  ee 
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(V '  V)  (vv) 


Now,  It  Is  evident  thet  if 
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yE  0  “  bE  0 
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\  <  1  FS  <  F 

^  -  1  then,  FS  -  F 
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FS>  FB 


By  definition, 


(B.12) 
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Vi  ’  bVt 
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Additionally,  for  all  cases  tabulated  in  Table  2, 

Fy  >  F* 


(B.16) 


Relation  of  Eq.  B.16  laplles  thet 


which  in  turn  laplles  thet 
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